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Principle  Investigator  Dr.  Karl  W.  Frese,  Jr. 

Cognizant  ONR  Scientific  Officer  Dr.  Harold  Guard 

Telephone  No.  4 1 5-859  322 1 

Project  Description 

The  goal  of  this  project  is  to  gain  an  understanding  of  the  factors  that  control  the  rate  of 
electrochemical  reduction  of  CO2  with  particular  emphasis  on  the  formation  of  methanol.  The 
approach  combines  theoretical  calculations  with  kinetic  measurements  on  various  metal  and 
semiconductor  surfaces.  The  theoretical  calculations  of  surface  thermodynamics  are  helpful  in 
identifying  allowed  elementary  steps  and  ruling  out  steps  that  are  thermodynamically  difficult. 
Mechanistic  determinations  are  aided  by  estimates  of  rate  constants  for  elementary  steps.  Therefore 
we  are  exploring  methods  of  calculating  rate  constants  for  elementary  steps  in  the  CO2  reduction 
pathway.  Laboratory  experiments  include  electrolysis  of  CO2  saturated  electrolytes  and  analysis  of 
product  distributions  as  a  function  of  electrode  identity,  surface  preparation,  electrode  potential, 
temperature,  and  pH. 

Significant.  Restills 

As  part  of  the  investigation  of  rate  constants  of  elementary  steps,  we  considered  the  addition  of  an 
electron  to  the  CO2  molecule  in  an  electrolyte.  The  activation  energy  for  the  uncatalyzed  process  is 
determined  by  the  reorganization  energy  of  the  molecule.  Two  methods  were  used  to  obtain  the 
reorganization  energy.  The  Marcus  equation  was  used  to  analyze  the  ab-initio  energy  surfaces  for 
the  gaseous  CO2/CO2'  couple.  This  gives  the  inner  reorganization  energy  of  2.5  eV,  a  large  value. 
The  electrostatic  contribution  was  found  to  be  0.6  eV  by  means  of  the  Marcus  equations.  The  total 
reorganization  energy  was  found  to  be  about  3. 1  eV.  The  major  reason  for  the  relatively  slow 
kinetics  has  been  found  to  be  the  bending  of  the  CO2  bond  angle  from  180  to  149°  in  the  transition 
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state.  Analysis  of  exchange  current  density  data  for  CO2  reduction  to  formate  (limited  by  electron 
addition  to  CO2)  yielded  3.2  eV  for  the  total  reorganization  energy. 

The  Sanderson  model  of  Polar  Covalence  was  proven  to  be  useful  for  the  calculation  of  the 
energetics  of  CO  adsorption  on  fifteen  metals.  The  heats  of  adsorption  of  various  stereochemical 
forms  including  linear,  bridged  and  dissociated  CO  were  successfully  calculated.  The  data  were 
compared  with  a  combination  of  HREELS,  LEED  and  thermal  desorption  data  from  the  literature. 
Average  deviation  of  the  calculated  results  from  the  experimental  values  was  about  2-3  kcal  mol'1. 
These  results  coupled  with  successful  calculations  of  binding  energies  of  H,  C  and  O  atoms  gives 
us  confidence  that  we  can  calculate  the  energetics  of  surface  reactions  involving  bond  with  transition 
metals  and  the  just  mentioned  atoms. 

We  also  developed  a  useful  method  for  predicting  the  bond  length  of  adsorbed  CO  molecules  on 
metal  surfaces.  The  method  is  based  on  quantum  mechanics  and  utilizes  the  partial  charges  on 
bonded  atoms  determined  by  the  Sanderson  method.  These  theoretical  data  were  needed  in  the 
above  calculations  of  CO  binding.  We  believe  they  are  accurate  to  0.005  A. 

The  rate  of  CO2  reduction  to  CH4,  CO  and  H2  was  studied  on  Cu  single  crystal  electrodes  in 
KHCO3  solution.  The  results  showed  that  the  (1 1 1)  plane  has  the  highest  methane  formation  rate. 
The  (100)  has  the  lowest,  with  (1 10)  intermediate.  This  is  the  first  evidence  that  crystal  plane  effects 
exist  for  CO2  reduction  on  metal,  as  opposed  to  semiconductor,  surfaces.  The  relative  rates 
correlated  with  the  heat  of  adsorption  of  CO.  The  CO  is  bound  weakest  on  (1 1 1)  and  strongest  on 
(100).  CO  has  been  shown  to  be  an  important  intermediate  in  the  CO2  to  CH4  pathway. 

Table  1  summarizes  results  on  the  rates  of  formation  of  methanol  on  various  electrode  surfaces  in 
aqueous  CO2  saturated  electrolytes.  Among  the  metals,  Mo  and  Pd  electrode  gave  the  highest  rates 
of  methanol  formation,  and  GaAs  (111)  the  highest  of  the  semiconductors  tested. 
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Table  1 

Materials  for  CO2  to  CH3OH  Conversion 


METAL 

ELECT. 

V(SCE) 

aRATE 

F 

GaAs(Aslll) 

Na2S04 

-1.2  to -1.4 

6  x  10-7 

100 

(100) 

0 

(HO) 

0 

InP(Plll) 

Na2S04 

-1.2  to -1.4 

3.7  x  10'7 

80 

Mo 

Na2S04 

-0.8 

5.0  x  10-7 

80 

Mo(cycled) 

Na2S04 

-0.8 

1.6  x  lO’6 

300 

Au 

KHCO3 

-1.2 

4.8  x  10-7 

38 

Ag  on  Cu 

KHCO3 

-1.35 

2.6  x  10-6 

63 

Re 

KHCO3 

-1.20 

1.9  x  10'6 

15 

10%  Ru/Pt 

KHCO3 

-1.20 

7.1  x  10'7 

9 

Cu(60°C) 

Na2S04 

-1.30 

2.0  x  10'6 

19 

Zn 

Na2S04 

-1.20 

9.7  x  10'7 

39 

Pd 

KHCO3 

-1.10 

1.8  x  10‘6 

59 

a.  mol  cirr^hr'l,  F  =  faradaic  efficiency 
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Emma  Fl  ans 

During  the  coming  year,  we  will  concentrate  on  three  areas.  We  will  continue  to  evaluate  the 
thermodynamics  of  all  conceivable  elementary  steps  in  the  CO2  reduction  pathway  on  various 
metal  surfaces.  The  goal  will  be  to  determine  the  most  likely  steps  and  suggest  how  they  may  vary 
on  different  metals. 

The  results  of  the  surface  thermodynamic  study  will  guide  us  towards  the  realization  of  an 
effective  catalyst  for  methanol  formation.  Electrodes  consisting  of  Cu/Zn  and  CuO/ZnO  will  be 
investigated  for  catalytic  activity  for  methanol  formation. 

We  will  apply  our  open  circuit  transient  technique  and  ac  impedance  methods  to  investigate 
overpotential  deposited  intermediates  on  electrode  surfaces  undergoing  CO2  reduction.  This 
information  is  necessary  for  the  elucidation  of  reaction  intermediates. 

We  have  received  a  renewal  contract,  extending  the  work  for  another  18  months.  During  this 
time  we  will  focus  on  the  role  of  oxides  and  composite  surfaces  in  CO2  reduction.  For  example, 
in  the  renewal  proposal,  we  suggested  a  strategy  for  ethylene  glycol  formation  on  Cu/Ag 
electrodes.  This  approach  utilizes  the  special  chemistry  at  Cu  to  form  ethylene  from  CO2  and 
oxygen  bearing  silver  to  catalyze  epoxidation  of  ethylene. 

Technical  Reports 

A  copy  of  the  published  version  of  technical  report  No  5  (  ref  1  above)  is  attached. 

Personnel 

Drs.  J.  J.  Kim  and  Mr.  Mike  Cinnibulk  also  contributed  to  the  ONR  program. 


9 


Reprinted  from  The  Journal  of  Physical  Chemistry,  1989,  93,  5911. 

Copynr'.t  ©  1989  by  the  American  Chemical  Society  and  reprinted  by  permission  of  the  copyright  owner. 


Calculation  of  Gibbs  Hydration  Energy  with  the  Ion-Dielectric  Sphere  Model 


Kart  W.  Frese,  Jr. 

Materials  Research  Laboratory,  SRI  International.  Menlo  Park,  California  94025 
( Received :  December  17,  1987;  In  Final  Form:  March  7,  19891 


The  Gibbs  energy  of  hydration  of  gaseous  ions  of  -  I  to  +4  charge  is  calculated  on  the  basis  of  the  exact  electrostatic  solution 
for  the  interaction  energy  of  an  ion  and  groups  of  dielectric  spheres.  Dispersion  energies  were  added.  The  average  error 
of  the  method  is  2%.  The  dielectric  constant  of  the  spheres  is  assumed  to  be  described  by  Booth’s  field-dependent  version 
of  Kirkwood's  theory  The  energetic  contributions  of  successive  layers  of  water  are  obtained  by  demanding  a  self-consistent 
dielectric  constant  at  each  shell  of  water.  The  results  confirm  that  78-90%  of  the  hydration  energy,  depending  on  the  ion. 
involves  water  in  the  inner  layer  The  accuracy  of  the  Born  expression  for  solvated  complex  ions.  M(H20)„*'\  is  in  doubt 
Inner-layer  hydration  numbers  ranging  from  4  to  12  were  obtained 


Introduction 

Model  calculation  of  the  Gibbs  energy  of  hydration  of  gaseous 
ions,  AGhyd,  has  a  long  history.1  Experimentally  derived  values2 3"6 
have  been  tabulated.  The  subject  has  been  reviewed  by  Case,4 
Friedman,  and  Krishnan,5  and  Conway.7  Survey  of  this  work 
shows  that  no  one  model  has  been  used  successfully  to  calculate 
AGhyd  for  ions  with  -2  to  +  4  charges.  Most  authors  have  limited 
their  comparisons  between  theory  and  experiment  to  univalent 
species.  For  example,  many  calculations8""  for  single-charge 
species  such  as  Na*  and  Cl"  ions  have  been  made  with  good  results. 
These  same  model  calculations,  some  of  which  utilize  the  Born 
model  outside  the  primary  hydration  layer  in  the  form  of  eq  la. 


AGbc 


(:eo)2 
2  (r,  +  rj 


Ua> 


are  usually  not  applied  to  higher  charged  species.  The  simple  Bom 
model  with  rw  =  0  is  known  to  fail  for  both  small  and  highly 
charged  ions.  Replacement  of  rw  with  empirical  values5  char¬ 
acteristic  of  univalent  anions  and  cations  reduces  the  errors  to 
about  10%  for  Fe2*  and  Fe3*.  Although  the  Born  model  is  ac¬ 
cepted  as  reliable  for  the  hydration  energy  of  the  complex  ion 
composed  of  a  univalent  ion  and  its  primary  hydration  layer,  it 
is  not  clear  that  such  a  model  would  be  accurate  for  +2,  +3,  or 
+4  complex  ions.  Furthermore,  the  electrostatic  model  treatments 
of  the  inner  layer  of  hydration  may  be  improved.  We  attempt 
to  show  below,  with  improved  models  of  the  inner-layer  hydration, 
that  eq  la  is  not  accurate  beyond  the  first  shell  of  water. 

A  number  of  studies  are  based  on  combination  of  eq  la  with 
the  dipole  and  quadrupole  moment  interaction  expressions  The 
model  of  Buckingham,8  for  example,  makes  use  of  ion-dipole,  eq 
lb.  lon-induced-dipole,  eq  lc.  and  ion-quadrupole  terms. 

(re0)Mw 

—  —TV  -  (lb) 

(r,  -I-  rw)‘ 
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Muirhead-Gould  and  Laidler,12  who  did  attempt  calculations  for 
-I  to  +4  ions,  also  included  questionable  lon-induced-dipole  terms. 
Generally,  ion-dipole  and  ion-induced-dipole  terms  are  calculated 
to  be  large  energies.  The  former  has  been  reported8 10  to  be  of 
the  order  of  the  heat  of  hydration.  Boucher13  in  1952  stated  that 
simple  expressions  such  as  eq  lb.c  are  not  accurate  at  ordinary 
bond  distances  of  2-3  A  in  the  inner  hydration  layer.  He  showed 
that  errors  of  50-100%  are  expected  Yet  these  expressions  were 
used  repeatedly  in  subsequent  years.  At  bond  distances,  the  r"1 
term  of  the  Born  model  is  an  overestimation,  and  the  r:  and  r~* 
terms  are  underestimations  of  the  true  energy  at  a  given  distance. 

In  addition,  the  quadrupole  moment  used  in  the  Buckingham 
model  is  about  an  order  of  magnitude  smaller414  than  the  measured 
value  for  water.  In  the  series  of  positive  ions  Li*.  Na*.  K*.  Rb*. 
and  Cs*.  the  magnitude  of  calculated10  ion-quadrupole  energy 
varies  from  53  to  26%  of  the  total  hydration  energy.  If  these  values 
are  revised  downward  an  order  of  magnitude,  the  agreement  with 
experiment  would  be  much  worsened.  As  an  additional  example 
of  the  failure  of  the  model  described  by  Buckingham,  consider 
the  Fe2*  ion.  The  calculated  heat  of  hydration  is  -676  kcal  mol"1 
compared  to  the  experimentally  derived  value2-5  of  -480  kcal  mol"1, 
an  error  of  about  40%.  Larger  errors  would  be  obtained  for  higher 
charges  and  for  smaller  ions  of  similar  charge  and  a  more  realistic 
quadrupole  moment. 

In  this  paper  we  present  an  improved  electrostatic  model  of 
ion-solvent  interactions  that  is  able  to  reproduce  the  hydration 
free  energies  of  ions  with  any  charge  type  -I  to  +  4  with  an 
average  error  of  2. 1  %.  The  ion-dipole,  ion-induced-dipolc.  and 
ion-quadrupole  terms  do  not  appear  explicitly  in  the  energy  ex¬ 
pression  The  contribution  of  each  shell  of  solvent  to  the  total 
hydration  energy  is  obtained.  The  equations  are  based  on  the  exact 
electrostatic  solution"  to  the  problem  of  the  interaction  energy 
of  a  charged  particle  with  a  neighboring  dielectric  sphere  The 
dielectric  constant, «,  of  the  spheres  depends  on  the  electric  field 
at  their  respective  positions.  The  use  of  the  dielectric  properties 
of  the  medium  is  analogous  to  the  dielectric  function  approach15 
used  by  phvsicists  in  the  theories  on  interaction  of  ions  and 
electrons  in  metals.  Quantum  mechanical  dispersion  energies16 


(12)  Muirhead-Gould,  J  S  ;  Laidler.  K  J  Chemical  Physics  of  Ionic 
Solutions,  Conwav.  B  E  .  Barradas.  R  G  .  Eds  .  Wilev  New  York.  1465, 
P  75 

(13)  Boucher.  C.  J  f  Theory  of  Electric  Polarization-.  Elsevier  Am- 
sierdam,  1952;  Chapter  5 

( 14)  Kistenmacher  H  .  Popkie,  H  .  Clementi.  E  J  Chem  Phvs.  1974.  61 . 
799 

(15)  Harrison.  W  A.  Electronic  Structure  and  the  Properties  of  Solids. 
W  H  Freeman  San  Francisco,  1980 

1 16)  Pi(7er.  K  S  Advances  in  Chemical  Physics:  Prigogme.  I  .  Ed  .  In- 
lerscience  New  York.  1959.  Vol  II,  p  69 
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are  added  to  the  electrostatic  results  to  give  the  chemical  part  of 
the  hydration  free  energy.  Bucher  and  Porter17  have  determined 
nonelectrostatic  short-range  interaction  energies  by  subtracting 
an  electrostatic  contribution  from  the  observed  Gibbs  energy  of 
hydration.  The  much  smaller  energy  associated  with  the  transfer 
of  the  ion  across  the  vacuum-solvent  interface5  was  not  considered 
in  the  present  work  and  poses  no  special  problem  because  ex¬ 
perimentally  derived  free  energies  of  hydration  do  not  include  this 
quantity. 

Our  motivation  for  this  work  is  related  to  an  important  question 
in  electrochemistry  concerning  the  role  of  solvent  molecules  held 
on  a  charged  electrode.  It  is  generally  assumed,  following  the 
detailed  calculation  of  Bockris  and  Hill18  on  Hg  electrodes,  that 
a  monolayer  of  solvent  covers  the  electrode.  This  means  that 
elementary  steps  in  certain  electrochemical  mechanisms  should 
be  considered  as  solvent  displacement  reactions19  and  therefore 
the  binding  energy  of  the  solvent  becomes  an  important  consid¬ 
eration.  Before  attempting  a  calculation  of  solvent  binding  on 
electrodes,  we  thought  it  necessary  to  seek  improved  methods  of 
calculation  of  the  energetics  of  electrostatic  binding  of  water  to 
both  ions  and,  by  analogy,  charged  electrode  surfaces. 


Theory 

Synthesis  of  an  Electrostatic  Model  for  SG^el).  The  BASIC 
computer  program  we  developed  for  the  electrostatic  contribution 
to  AGhy*  is  based  upon  an  equation  from  Bottcher'3  that  quantifies 
the  interaction  energy  between  an  ion  of  charge  q  and  a  neighbor 
spherical  region  or  regions  of  dielectric  constant,  t.  The  term  W 
below  represents  the  electrostatic  free  energy  change13  for  the 
process  in  which  the  ion  and  sphere  are  brought  together  from 
infinite  separation.  In  our  calculations  the  neighboring  spherical 
regions  correspond  to  water  molecules  spatially  ordered  by  in¬ 
teraction  with  the  ion  and  by  hydrogen  bonding.  Thus  the  model 
is  structural  rather  than  continual.  The  relationship  given  by 
Bottcher  for  the  interaction  energy  between  an  ion  of  charge,  q, 
at  a  distance,  s.  from  a  single  dielectric  sphere  of  radius  o  is 


-W=  — - - -  - 

2  /-o</  +  /  +  1  5 


2/+1 
21+  2 


(2) 


Here  a  =  rw  =  1 .38  A,  the  radius  of  the  water  molecule  as  deduced 
from  X-ray  and  electron  diffraction  data20  on  ice.  The  quantity 
s  =  r,  +  rw,  where  rx  is  the  ionic  radius.  In  general  ^  =  R(K)  as 
discussed  below.  Equation  2  follows  from  a  solution  of  Laplace's 
equation  under  appropriate  boundary  conditions.13  Note  that  the 
usual  properties  ascribed  to  solvent  molecules,  such  as  dipole  or 
quadrupole  moments,  do  not  appear,  but  rather  the  general  di¬ 
electric  constant  function.  Clearly,  «  will  depend  on  the  dipole 
moment,  m.  and  other  molecular  properties  such  as  the  Kirkwood 
correlation  factor,21  g,  which  is  closely  related  to  the  solvent 
structure.  The  dielectric  constant  in  eq  2  also  depends  on  the 
electric  field  strength. 

The  second  building  block  in  the  model,  eq  3,  due  to  Booth21-22 
provides  the  necessary  expression  for  «  as  a  function  of  electric 
field  strength  and  thus  distance  from  the  ion: 


(  =  n2  + 


28iV07r(n2  +  2)m, 


3(73) 


coth  ( x ) 


(3) 


where 


x 


731/2  fjjn2  +  2 )E 

6  kr 
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and  the  term  in  brackets  is  the  Langevin  function.  In  eq  3  ,V0 
is  the  particle  density  of  water  (cm-3);  fiy,  the  dipole  moment 
according  to  Booth,  is  2. 1  X  I  O'18  esu-cm;  n  is  the  optical  refractive 
index;  and  £  is  the  electric  field  strength  in  esu  Equation  3,  which 
is  based  on  Kirkwood’s  theory,  contains  Booth’s  approximate 
correlation  factor  derived  from  the  hydrogen-bonded  structure 
of  water.  Equation  3  has  been  used  in  several  previous  studies 
of  hydration.4 

The  next  step  in  assembling  the  theory  was  to  generalize  eq 
2  and  3  so  that  e  and  -  W  could  be  calculated  at  the  positions  of 
the  water  molecules  at  various  distances  from  the  central  ion.  The 
field  around  an  ion  is  a  function  of  distance;  therefore,  we  make 
the  calculation  of  a  “differential  hydration  energy"  and  then  sum 
the  results  at  each  distance.  This  part  of  the  calculation  recognizes 
the  discrete  structure  of  liquid  water  (a  hydrogen-bonded  network) 
and  is  consistent  with  the  familiar  notion  of  various  layers,  shells, 
or  spheres  of  hydration,  i.e.,  primary,  secondary,  and  so  forth 
Using  this  method  we  calculated  the  contribution  to  AGhyd(el)  of 
various  shells  of  water  molecules  as  a  function  of  distance.  Then 
the  contributions  from  each  shell  at  a  certain  distance  are  summed 
to  give  the  final  result  for  the  electrostatic  contribution. 

The  input  to  the  computer  program  consisted  of  2,  the  integral 
charge  on  the  ion;  N,  the  number  of  water  molecules  in  the  primary 
hydration  layer;  and  r„  the  crystal  radius  of  the  ion  as  given  by 
Pauling.23  Marcus'  has  concluded  that  the  Pauling  radius  is  a 
good  measure  of  the  size  of  solvated  ion.  He  showed  that  a  plot 
of  crystal  radius  versus  ion-water  distances  derived  by  a  variety 
of  methods  incuding  Monte  Carlo,  molecular  dynamics.  X-ray, 
and  neutron  diffraction  gave  a  linear  relation  with  unit  slope  and 
intercept  of  1.39  A  for  the  radius  of  the  water  molecule. 

The  formula  (in  A)  for  the  distance  between  the  central  ion 
and  the  Kth  dielectric  sphere  is 

R(K)  =  (r,  +  1.38)  +  2.76K  (4) 

where  K  was  varied  from  0  to  9.  K  =  0  corresponds  to  the  first 
or  inner  hydration  shell.  The  dielectric  constant  is  given  by 
j\E{  *?(£)]),  where / stands  for  Booth’s  eq  3.  The  computational 
task  was  to  find  {[/?(£)]  at  each  of  the  Afth  shells  of  water 
molecules.  This  was  accomplished  by  first  assuming  t  =  I  and 
then  calculating  E(K)  as  21I\R(K)  t(K)]  This  first  estimate  of 
the  field  gave  the  second  estimate  of  e  by  means  of  eq  3  and  thus 
a  second  estimate  of  the  Held  at  R(K).  This  iterative  process  was 
performed  until  self-consistency  was  achieved  at  each  R(K).  The 
criterion  for  terminating  the  calculation  loop  was  that  two  con¬ 
secutive  t(AO  differed  by  less  than  0.1%.  The  resulting  value  <(A1 
was  used  in  the  first  10  terms  in  the  summation  over  /  in  eq  2  to 
give  the  contribution  to  AGhyd(el)  at  that  distance.  The  number 
of  ions  in  the  primary  hydration  shell  (K  =  0)  is  .V  as  given  in 
Table  VII.  At  each  distance  beyond  the  first  hydration  shell.  K 
>  0,  the  number  of  neighbors  was  assumed  to  be  given  bv  NIK) 
=  INK. 

This  algorithm  follows  from  the  assumption  that  each  water 
molecule  in  the  primary  shell  can  bind  two  additional  water 
molecules  by  means  of  hydrogen  bonding,  thus  forming  the  second 
shell.  The  molecular  binding  sites  are  either  the  two  lone  pairs 
of  electrons  on  oxygen  or  the  pair  of  hydrogen  atoms.  The  number 
of  neighbors  in  the  third  shell  is  similarly  composed  of  the  sum 
of  the  pairs  of  water  molecules  bound  to  each  molecule  in  the 
second  shell  The  exact  prescription  for  calculating  the  number 
of  neighbors  beyond  the  second  shell  is  relatively  unimportant 
because  of  the  small  energy  contribution  from  far  shells.  In  the 
calculations,  the  first  10  discrete  distances  of  water  spheres  from 
the  central  ion  were  included  to  estimate  the  energetic  contribution 
of  the  roughly  first  30  A  of  water  around  the  infinitely  dilute  ion 
This  means  that  10  spheres  of  hydration  were  calculated  Fol¬ 
lowing  Buckingham,810  the  hard-sphere  model  was  assumed  for 
the  repulsive  energy. 

Quantum  Mechanical  Dispersion  Energies.  The  chemical  free 
energy  of  hydration.  A Ghyd.  was  calculated  by  adding  the  dispersion 


(23)  Pauling.  1.  The  Mature  of  the  Chemical  Bond.  3rd  ed  ,  Cornell 
University  Press  Ithaca  NY.  I960;  Chapter  I  3 
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TABLE  I:  Results  for  Ns* 


TABLE  V:  Results  for  Fe3* 


R(K),  A 

E.  esu/cm 

< 

(cos  8) 

-W(K).  eV 

R(K).  A 

E.  esu/cm 

< 

(cos  8) 

-W[K).  eV 

2.33 

38200 

23.2 

0.903 

3.51 2 

2.02 

1486000 

2.38 

0997 

38.887 

;.09 

2266 

81.8 

0.200 

0.235 

4.78 

11170 

56.5 

0.675 

4,054 

27.17 

77 

83.8 

0.007 

0.001 

7.54 

3168 

80  0 

0  274 

0945 

3.858“ 

26.86 

238 

83.8 

0021 

0018 

-AG^fel). 
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TABLE  II:  Results  for  Rb* 


R(K).  A 

£,  esu/cm 

t 

(cos  9) 

-W(K),  eV 

2.86 

9669 

60.7 

0.630 

2.462 

5.62 

1843 

82.4 

0.164 

0.273 

27.7 

75 

83.8 

0.007 

0.002 

2.891“ 

•-AGhyd(el). 

TABLE  III:  Results  for  F" 

R(K).  A 

£,  esu/cm 

(cos  8) 

eV 

2.74 

11530 

55.5 

0.684 

3.419 

5.50 

1928 

82.4 

0.171 

0.341 

27.58 

75 

83.8 

0.006 

0.002 
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TABLE  IV:  Results  for  Ft1* 


/?(*),  A 

£,  esu/cm 

e 

(cos  8) 

-W(K).  eV 

2.14 

677000 

3.10 

0.994 

16.086 

4.9 

5409 

73.9 

0.431 

1.644 

7.66 

1990 

82.3 

0.177 

0.394 

26.98 

157 

83.8 

0.014 

0.008 

18.452“ 

*-AGw(e  I). 

energies  to  the  classical  electrostatic  energy,  AGhyd(el),  outlined 
above.  We  followed  the  discussion  of  Pitzer16  and  employed  the 
Kirkwood-Muller  dispersion  energy  formula,  eq  5,  for  a  pair  of 
particles  A  and  B. 


ft  =  _ 
disp 


hm<faKa B|  aK  aB 
No#  I  Xa  Xb  , 


The  dispersion  energy  depends  on  the  sixth  power  of  the  separation, 
the  electronic  polarizability,  a,  and  diamagnetic  susceptibility, 
X,  of  the  particles.  The  appropriate  physical  constants  have  been 
tabulated16  for  each  inert  gas  configuration.  We  assumed  that 
the  ion-water  dispersion  interaction  could  be  represented  by  the 
appropriate  inert  gas-neon  interaction.  The  neon  properties  were 
chosen  to  represent  the  water  molecule  because  most  of  the 
electron  density  in  the  molecule  resides  on  the  oxygen  atom. 
Inserting  the  physical  constants  and  using  Ne  properties  to  rep¬ 
resent  the  water  molecule,  we  obtain  eq  6,  expressed  in  kcal  mol'1 


I'd  nr 


46.92JVaB 


“B 

0.0555  +  — 


where  B  represents  any  solvated  ion.  In  eq  6,  distance  is  expressed 
in  angstroms,  polarizability  is  expressed  in  I  O'24  cm3,  susceptibility 
is  expressed  in  10"6  cm3,  and  N  is  the  number  of  water  molecules 
in  the  primary  hydration  sphere.  The  ratio  of  the  magnetic 
y  Ability  to  ihe  polarizability  for  the  ions  was  calculated  from 
f  propriate  inert  gas  values.16  If  the  ion  did  not  have  an  inert 
►  .lfiguration,  interpolated  or  extrapolated  (in  the  case  of 
.-like  ions)  values  were  employed.  For  example,  the  properties 
of  Ft'  'ere  interpolated  between  the  values  for  Kr  and  Ar 
acc  _  to  the  number  of  electrons. 

Tn  auling  crystal  radii  for  certain  transition-metal  ions  were 
■’ssumed  to  reflect  the  crystal  field  stabilization  energy  (crystal 
bond  shortening),  and  therefore  no  crystal  field  energies  were 


TABLE  VI:  Results  for  Ce4* 

R(K).  A  £,  esu/cm 
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(cos  9) 
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VII:  Comparison  of  Calculated  and  Observed  Chemical  Free 
of  Hydration  of  Caseous  Ions  at  25  °C  (kcal  mol1)* 


ion 

rc.  A 

N 

AGdup 

-AGhyd(calcd) 

-agm‘ 

-agm' 

H* 

0.3924 

4 

0 

(260.1) 

260  5 

259.2 

Li* 

0.78 

4 

3.9 

119  6 

123.5 

120  8 

Na* 

0.95 

4 

4.3 

93.2 

98.3 

9'7,0 

Ag* 

1.26 

8 

2.0 

108 

1 14.5 

113 

K* 

1.33 

6 

5.5 

77.0 

80  8 

79  3 

Rb* 

1.46 

7 

9.1 

75.8 

76.6 

'4  2 

Cs* 

1.69 

8 

12 

69.2 

71.0 

66.5 

F" 

1.36 

8 

3.2 

94.2 

103.8 

88.2 

OH' 

1.40' 

8 

2.9 

88.6 

90  6 

CL 

1.81 

12 

6.2 

79  9 

75.8 

74.8 

Br‘ 

1.95 

12 

7.1 

69.4 

72.5 

67.9 

I- 

2.16 

12 

7.7 

57.1 
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59.0 

Mg2* 

0.64 

5 

12 

454 

455 

453 

Cu2* 

0.72 

6 

38 

497 

499 

495 

Zn2* 

0.74 

6 

44 

486 

485 

482 

Fe2* 

0.76 

6 

42 

467 

456 

450 

Ca2* 

0.99 

7 

21 

365 

381 

378 

Cd2* 

0.97 

8 

39 

443 

431 

428 

Hg2* 

1.10 

9 

60 

443 

434 

Pb2* 

1.20 

8 

53 

355 

358 

355 

Ba2* 

1.35 

9 

27 

318 

315 

313 

A)3* 

0.50 

4 

15 

997 

1103 

1099 

Fe3* 

0.64 

6 

49 

1079 

1036 

1032 

Cr3* 

0.69 

6 

53 

1077 

1037 

1036 

Eu3* 

1 .03 

9 

59 

815 

843 

851 

Ce3* 

1.11 

10 

53 

805 

804 

808 

!  a3* 

1.15 

10 

48 

761 

791 

777 

Ce“* 

1.01 

10 

68 

1448 

1471' 

“r,  is  the  Pauling  crystal  radius.  ,V  is  the  primary  hydration  num¬ 
ber.  -AGhyd(calcd)  is  the  sum  of  -AGhyd(el)  and  -AGd„p.  Chemical 
free  energies  means  that  the  small  energy  change  associated  with  the 
potential  drop  across  the  vacuum-solvent  interface  is  not  included  in 
the  hydration  free  energies  in  this  paper.  ‘Reference  4.  experimental 
values,  'Reference  1.  experimental  values.  4  Empirical  radius  to  give 
to  -260.1  kcal  mol-1.  'Thermochemical  radius,  ref  24.  >  Data  from  ref 
24  and  43. 

added  to  the  calculated  hydration  energies. 

Results  and  Discussion 

Typical  results  for  the  dielectric  constant  and  hydration  energy 
as  a  function  of  field  and  distance  are  given  in  Tables  I— V I  for 
various  ions.  The  calculated  average  cosine  of  the  angle  between 
the  field  and  the  water  dipole  is  also  given.  These  values  show 
the  range  of  influence  of  the  electric  field  of  the  ion.  For  singly 
charged  ions,  (cos  8)  falls  to  0.1  at  8-10  A  and  1 0-1  5  A  for  +3 
or  +4  species.  The  free  energy  value  at  the  bottom  of  the  last 
row  represents  the  summation  of -IT  =  AGhyd(el)  at  each  distance. 

The  Gibbs  free  energies  of  hydration,  given  in  the  fourth  column 
of  Table  VII,  were  obtained  by  adding  the  dispersion  energies 
(column  7)  to  AGhyd(el).  Columns  5  and  6  contain  experimentally 
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derived  values.  In  general  the  overall  agreement  is  very  good. 
Note  that  the  effective  radius  of  the  proton,  0.392  A,  is  the  value 
required  to  obtain  agreement  with  experiment,  assuming  N  =  4. 
Note  that  Bucher  and  Porter17  found  0.35  A  by  a  different  type 
of  electrostatic  calculation.  The  OH'  radius,  1.40  A,  is  the 
thermochemical  radius24  used  in  lattice  energy  calculations.  We 
used  Pauling  crystal  radii  for  all  other  ions  with  the  exception  of 
Li+  for  which  the  Goldschmidt  value,1'25  0.78  A,  was  used.  The 
inapplicability  of  the  standard  Pauling  radius,  0.60  A,  can  be 
explained  by  analogy  to  the  situation  in  Li+  salts.  Due  to  the  small 
size  of  the  Li+  ion,  the  anion-anion  contact  repulsions  play  a 
critical  role  in  keeping  the  lithium-anion  distance1  greater  than 
the  sum  of  the  standard  Pauling  radii.  By  analogy,  intermolecular 
water  molecule  repulsions  are  operational  in  the  primary  solvation 
shell  causing  a  similar  effect.  We  are  not  aware  of  any  accurate 
value  for  the  free  energy  of  hydration  of  S2';  therefore  sulfide  ion 
was  not  included. 

The  average  percent  deviations  between  the  calculated  and 
experimentally  derived  results  for  each  charge  are  the  following: 
-I,  3.5%;  +1,  2.5%;  +2,  1.7%;  +3,  1.5%.  In  the  single  +4  case 
of  Ce4+,  the  deviation  was  1.5%.  The  calculated  results  are 
proportional  to  N ,  the  primary  hydration  number.  This  means 
that  increasing  or  decreasing  the  N  value  by  1  would  lead  to  a 
relative  deviation  in  the  calculated  AGhyd  of  (100 /N)%.  The 
relative  error  is  25%  for  N  =  4  and  1 2%  for  N  =  1 2.  This  deviation 
is  larger  than  the  experimental  uncertainty  in  A<7hyd,  and  therefore 
the  integral  hydration  numbers  are  unique  within  the  model. 

Primary  Hydration  Numbers.  A  large  body  of  data  on  hy¬ 
dration  numbers  is  presented  in  the  review  by  Hinton  and  Amis.25 
A  well-known  difficulty  is  that  of  finding  agreement  among  the 
various  experimental  methods  for  deriving  hydration  numbers. 
We  decided  against  detailed  comparison  of  our  N  values  to  those 
in  the  review25  because  of  the  ambiguity  in  the  literature.  This 
should  not  be  taken  to  mean  that  our  A  values  are  unreasonable. 
Certainly,  we  could  justify  most,  if  not  all,  of  our  N  values  by 
quoting  results  obtained  by  several  different  methods.  But,  it  is 
doubtful  that  a  single  method  would  yield  results  consistent  with 
ours  for  all  ions.  Indeed  no  single  method  has  been  applied  to 
all  the  ions  in  Table  VII. 

Some  of  the  reasons  for  the  difficulties  one  encounters  are  the 
following.  A  number  of  methods  depend  on  assignment  of  a 
standard  hydration  number  to  a  reference  ion.  We  found  it 
difficult  to  find  a  set  of  results  that  gave  both  a  reference  value 
and  the  sought-for  value  consistent  with  our  results.  Different 
methods  do  not  even  agree  on  whether  the  hydration  number 
should  increase  or  decrease  with  the  size  of  the  ion.  For  example, 
Conway7  has  noted  that  primary  hydration  numbers  determined 
by  nuclear  magnetic  resonance26-27  increase  with  ionic  radius 
whereas  values  derived  from  other  methods  including  mobility25 
indicate  smaller  hydration  numbers  for  larger  ions.  To  compare 
our  N  values,  one  must  select  a  method  that  measures  the  primary 
hydration  sphere.  This  generally  means  NMR  and  X-ray  dif¬ 
fraction28  methods  in  concentrated  electrolytes.  Experimental29 
and  theoretical30  evidence  indicates  that  apparent  hydration 
numbers  decrease  with  increasing  salt  concentration.  Our  model 
and  most  model  calculations  to  date  apply  to  infinite  dilution. 
Further,  there  is  evidence51  that,  in  concentrated  solutions,  the 
primary  hydration  number  depends  on  the  identity  of  the  coun¬ 
terion. 

The  clear  trend  in  our  results  is  that  primary  hydration  numbers 
increase  with  the  size  of  the  ion.  Studies  that  have  indicated  larger 
primary  hydration  numbers  for  larger  ions  include  NMR,26-27 
molecular  dynamics30"33  and  ab  initio  quantum14  calculations,  and 
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ION  RADIUS  |A| 

Figure  1.  Observed  ligancy  in  oxides  (solid  line)  and  inner-sphere  hy¬ 
dration  numbers  found  in  this  work  versus  Pauling  crystal  radius. 


X-ray  diffraction.28  By  analogy  to  ligancy  trends  in  crystalline 
solids,  we  might  expect  larger  primary  hydration  numbers  for 
larger  ions.  The  actual  number  of  ligands  in  the  solid  state  is 
mainly  determined  by  radius  ratio  effects.23 

We  looked  for  a  correlation  between  our  N  values  and  the 
number  of  nearest  neighbors  expected  on  the  basis  of  crystal 
packing.  The  analogy  between  frozen  water  in  the  primary  hy¬ 
dration  shell  and  solid-state  water,  ice,  has  been  made  before. 
Oxides  were  chosen  for  comparison,  and  Pauling23  has  tabulated 
the  expected  and  observed  ligancies  as  a  function  of  ionic  radius. 
The  observed  values  versus  the  known  ligancy  are  represented  by 
the  solid  line  in  Figure  1.  The  points  represent  the  values  of  N 
(Table  VII)  found  in  the  present  calculations.  It  is  seen  that  the 
observed  ligancies  in  the  solid  state  and  in  the  highly  polarized 
liquid  state  in  the  inner  hydration  layer  correlate  well,  especially 
for  univalent  ions.  Our  N  values  are  integers,  again  in  analogy 
to  solid-state  ligancy  and  in  contrast  to  many  fractional  experi¬ 
mental  determinations.25  There  is  little  controversy  concerning 
the  primary  hydration  number  of  four25  for  H+,  Li+,  and  Na*. 
The  N  =  4  for  Al3+  does  not  agree  with  the  commonly  accepted 
value34  of  6.  The  transition  elements  Fe2+,  Fe3+,  Zn2+,  and  Cu2+ 
were  required  to  be  hexaaquo  complexes35  to  obtain  agreement 
with  experimental  hydration  energies.  We  found  nine  inner-layer 
water  molecules  for  Eu3+;  a  value  of  eight  or  nine  is  supported 
by  spectroscopic36-37  evidence.  The  value  of  ten  for  the  three  other 
lanthanides38  in  Table  VII  is  very  reasonable  given  the  rather 
confident  value  of  nine  for  Eu3+. 

The  most  unconventional  result  is  N  -  1 2  for  the  halogens  Br‘, 
Cl',  and  I".  Such  large  numbers  are  consistent  with  the  icosahedral 
arrangement,23-39  which  represents  the  largest  number  of  spherical 
particles  that  can  be  packed  around  an  isodimensional  central 
sphere.  As  a  group,  these  halogen  ion  values  fall  very  close  to 
the  solid-state-ligancy  line  in  Figure  1,  yet  hydration  numbers 
between  4  and  0  have  been  reported.25  The  following  values  have 
been  determined  for  P,  Cl",  Br",  and  I",  by  the  NMR  method26 
and  based  on  N  -  6  for  K+;  N  -  9.9,  13.2,  16.2,  and  21.8. 
respectively.  Unfortunately,  N  =  1  was  found  for  Li+  in  the  same 
investigation.  Molecular  dynamics  calculations31  give  N  =  8-9 
for  Cl"  in  2.2  M  CsCl.  X-ray  diffraction  values28  for  the  inner 
layer  range  from  6  to  9  and  may  depend  on  concentration.  These 
results  support  our  large  inner-sphere  solvation  numbers  and  at 
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Calculation  of  Gibbs  Hydration  Energy 


TABLE  VIII:  Comparison  of  Inner-Sphere  Free  Energies  of 
Hydration  of  Gascons  loss  at  25  °C  (heal  mot*1) 


ion 

N 

-AG-wt  inner)* 

-(A H  -  7AS)‘ 

*A<W 

-A(/3C 

H* 

4 

248 

231 

260.1 

52 

Li* 

4 

111 

92 

119.6 

46 

Na* 

4 

85 

78 

93.2 

44 

K* 

6 

68 

73 

77.0 

40 

F* 

8 

82 

86 

94.2 

40 

OH* 

8 

76 

86 

88.6 

39 

•Results  of  this  work,  including  dispersion  energy;  see  Table  I  and 

VII  for  examples.  ‘Enthalpy  from  either  Hartree-Fock  (H-F)  calcu¬ 
lations  (ref  14)  or  experimental  data  on  gas-phase  cluster  reactions  (ref 
41  and  42).  Total  AS  of  hydration  (ref  5).  See  text  for  further  ex¬ 
planation.  'Total  calculated  Gibbs  energy  of  hydration,  Table  VII. 
< Gibbs  energy  of  hydration  for  outer  sphere  calculated  using  eq  la. 
Ionic  radii  from  Table  VII;  rw  *  1 .38  A. 

the  same  time  point  out  the  state  of  confusion  about  primary 
hydration  numbers. 

Inner -  and  Outer-Layer  Gibbs  Energy  of  Hydration.  We  define 
the  inner  (primary)  hydration  layer  to  include  the  N  water 
molecules  assumed  to  be  in  closest  contact  with  the  ion.  The  outer 
hydration  sphere  includes  water  molecules  beyond  the  primary 
sphere.  These  definitions  are  similar  to  those  employed  in  the¬ 
ories40  of  electron-transfer  rates.  It  is  frequently  assumed810  that 
the  outer  hydration  free  energy  of  the  ion-primary  water  complex 
is  described  by  the  Born  model,  eq  la.  The  radius  of  this  ion  is 
given  by  the  crystal  radius  of  the  ion  plus  2.76  A,  the  diameter 
of  a  water  molecule.  Using  the  results  in  Tables  I— VI,  we  can 
determine,  according  to  our  calculations,  how  much  of  AGhyd(el) 
is  attributed  to  ion-water  interactions  outside  the  primary  sphere. 
For  example,  the  calculated  electrostatic  interaction  energy  of 
Na(H2Oj4+  with  the  remaining  solvent  is  3.86  -  3.51  =  0.33  eV. 
Therefore  90%  of  the  hydration  energy  of  the  Na+  resides  in  the 
inner  layer.  Similar  figures  for  P,  Fe3+,  and  Ce4+  are  86,  87, 
87,  and  78%,  respectively.  The  smaller  percentage  for  Ce4+  results 
from  the  larger  contribution,  14%,  from  the  secondary  hydration 
sphere.  Application  of  the  Bom  model  to  the  Na(H20)4+  complex 
ion  gives  AG^outer)  =  1 .83  eV,  a  value  5  times  larger  than  found 
in  the  present  calculations.  Such  a  discrepancy  clearly  signifies 
a  major  fault  in  one  of  the  models. 

Our  calculations  show  that  78-90%  of  the  electrostatic  energy 
arises  in  the  primary  hydration  sphere  and  that  application  of  the 
3orn  model  gives  a  large  overestimation  of  the  hydration  energy 
of  the  ion-primary  water  complex.  This  result  is  not  expected 
from  previous  electrostatic  models.4-5-8*10 

The  conclusion  that  most  of  the  hydration  energy  'esidues  in 
the  inner  layer  follows  from  both  quantum  mechanical  calculations 
and  experimental  energetics  of  ion-water  cluster  reactions.  Table 

VIII  shows  data  that  support  our  results.  The  third  column  gives 
our  calculated  Gibbs  energy  of  hydration  for  the  inner-layer 
complex  ion  including  the  dispersion  energies.  The  appropriate 
values  for  Na+,  for  example,  may  be  obtained  from  Tables  I  and 
VII.  As  noted  above,  these  numbers  represent  80-90%  of  the  total 
energy.  The  fourth  column  in  Table  VIII  contains  an  approxi¬ 
mation  to  the  same  Gibbs  energy  based  on  the  ab  initio  Har¬ 
tree-Fock  (H-F)  calculations14  or  the  experimental  cluster  en¬ 
thalpies41-42  as  follows.  The  experimental  cluster  enthalpy  for 
H(H20)4+  was  used.  For  Li(H20)4+,  Na(H20)4+,  and  K(H20)6+, 
the  H-F  results  were  used;  these  were  in  close  agreement14  with 
the  experimental  cluster  enthalpies.  In  the  case  of  F(H20)s*,  the 
H-F  results  for  F(H20)5*  were  supplemented  for  the  sixth  to 
eighth  water  molecules  by  extrapolation  of  the  H-F  energy  vs  N 
plots  (see  Figure  13,  ref  14).  The  result  was  an  additional  sta¬ 
bilization  of  14  kcal  mol*1  for  the  water  molecules  6-8.  The 
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enthalpy  data  for  0H(H20)8*  were  obtained  by  also  adding  14 
kcal  mol"1  to  the  experimental  cluster  enthalpy41  of  OH(H20)5*. 
The  similar  treatment  of  solvated  P  and  OH"  is  justified  because 
they  are  isoelectronic. 

In  all  cases,  the  total  observed  entropy5  of  hydration  of  the 
gaseous  ion  was  added  to  the  H-F  or  cluster  enthalpy  on  the 
assumption  that  all  the  entropy  change  is  attributable  to  inner 
layer.  This  procedure  makes  the  Gibbs  energies  in  column  4 
uncertain  by  several  kcal  mol"1,  because  most,  but  not  all,  the 
entropy  change  is  attributable  to  the  inner  layer. 

The  last  column  in  Table  VIII  gives  the  Born  contribution  to 
the  Gibbs  energy  of  hydration,  eq  la,  for  the  inner-layer  complex 
of  the  ion  and  N  water  molecules.  The  agreement  between 
columns  3  and  4  of  Table  VIII  gives  strong  support  to  the  cal¬ 
culation  method  based  on  our  eq  2  and  3  above  and  the  fact  that 
inner-layer  solvation  dominates  by  far.  The  difference,  column 
5  minus  column  3,  is  seen  to  be  small  compared  to  the  Born 
charging  term.  1  his  supports  the  conclusion  that  the  energetic 
contributions  to  the  hydration  energy  beyond  the  inner  layer  are 
much  smaller  than  eq  la  would  lead  one  to  believe. 

As  stated  by  Bottcher,13  the  simple  expressions,  like  eq  lb,c, 
used  in  the  description  of  the  inner  hydration  layer,  underestimate 
the  interaction  energy  by  50-100%.  These  equations  are  based 
on  conducting,  point  ion  models,  which  are  not  valid  at  distances 
such  that  the  radius  of  the  particles  is  of  the  order  of  their  sep¬ 
aration.  If  this  is  correct,  it  would  explain  how  the  Born  model, 
eq  la,  in  combination  with  eq  lb,c  for  solvation  of  the  ion-primary 
water  complex  was  perceived  to  give  accurate  results.  The  un¬ 
derestimation  of  the  Gibbs  energy  of  hydration  in  the  inner  layer 
is  canceled  by  the  overestimation  in  the  outer  region  (outside  the 
first  shell  of  water),  and  the  sum  of  the  energy  in  the  two  regions 
is  approximately  correct,  at  least  for  univalent  ions. 

These  observations  may  have  some  bearing  on  the  activation 
process  for  homogeneous  and  heterogeneous  electron  transfers  to 
solvated  ions.  There  has  been  disagreement  concerning  the  relative 
importance  of  inner-  and  outer-sphere  (-layer)  contributions  to 
the  activation  energy.  According  to  Bockris  and  Kahn,40  the 
outer-sphere  activation  would  require  simultaneous  participation 
of  a  large  number  of  solvent  molecules,  each  having  a  small 
contribution  to  the  activation  energy.  It  is  stated  that  this  si¬ 
multaneous  activation  is  highly  improbable.  The  reorganization 
of  the  equilibrium  solvation  coordinates  is  generally  held  to  be 
the  rate-determining  step  in  the  electron  transfer.  Our  results 
suggest  that  the  solvation  energy  resides  mainly  in  the  inner 
hydration  layer,  at  least  in  the  ideal  solution  state.  This  may  mean 
that  the  reorganization  process  involves  a  relatively  few  solvent 
molecules,  the  order  of  N  in  Table  VII.  The  contribution  of  each 
ion-solvent  “bond"  in  an  exchange  reaction42  would  then  be  the 
order  of  A/4 N,  where  A  is  the  observed  reorganization  energy. 
The  change  in  inner-layer  hydration  with  concentration  and 
counterion  mentioned  above  may  be  seen  in  the  kinetics  of  electron 
transfer. 

Summary.  The  Gibbs  hydration  energy  has  been  calculated 
for  30  ions  with  -1  to  +4  charge  with  typical  deviations  of  2% 
from  the  experimentally  derived  values.  The  inapplicable  Born- 
type  equations  and  other  simple  electrostatic  functions  were  re¬ 
placed  by  a  more  accurate  series  solution  to  the  electrostatic 
problem  of  the  interaction  of  a  ion  and  a  dielectric  sphere.  The 
classical  electrostatic  part  of  the  calculation  involves  only  the 
dielectric  constant  and  discrete  structure  information.  Quantum 
mechanical  dispersion  energies  were  added.  Primary  hydration 
numbers  were  chosen  to  give  agreement  between  calculated  and 
experimentally  derived  hydration  energies.  Small  integer  values 
ranging  from  4  to  1 2  were  obtained.  These  values,  especially  the 
smaller  ones,  are  unique  to  the  model.  Pauling  ionic  radii  were 
employed,  except  for  the  Li+,  where  ligand  repulsions  are  probably 
operational  and  the  empirical  Goldschmidt  value  gives  much  better 
results.  The  resulting  primary  hydration  numbers  were  correlated 
with  the  ionic  radius  as  may  be  expected  from  known  ligancy-size 
trends  in  solid-state  chemistry.  The  calculations  suggest  that  not 
only  does  the  Born  model  not  generally  describe  the  primary 
solvation  shell  but  also  it  does  not  hold  for  the  solvation  energetics 
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of  the  complex  ion  formed  by  the  central  ion  and  the  primary 
hydration  shell. 

Acknowledgment.  This  work  was  supported  by  the  Chemistry 
Division  of  the  Office  of  Naval  Research. 

Registry  No.  Fi*.  12586-59-3;  Li*.  17341-24-1;  Na*.  17341-25-2; 


Ag*,  14701-21-4;  K*.  24203-36-9;  Rb*,  22537-38-8;  Cs*.  18459-37-5; 
F\  16984-48-8;  OH',  14280-30-9;  CL.  16887-00-6;  Br\  24959-67-9;  I", 
20461-54-5;  Mg2*.  22537-22-0;  Cu2*.  15158-11-9;  Zn2*,  23713-49-7; 
Fe2*,  15438-31-0;  Ca2*,  14127-61-8;  Cd2*.  22537-48-0;  Hg2*.  14302- 
87-5;  Pb2-,  14280-50-3;  Ba2*,  22541-12-4;  Al3*.  22537-23-1;  Fe3*. 
20074-52-6;  Cr3*,  16065-83-1;  Eu3*,  22541-18-0;  Ce3*,  18923-26-7; 
La1*.  16096-89-2;  Ce4*,  16065-90-0 


